Ginzburg-Landau theory of the cluster glass phase 
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On the basis of a recent field theory for site-disordered 
spin glasses a Ginzburg-Landau free energy is proposed to de- 
scribe the low temperatures glassy phase(s) of site-disordered 
magnets. The prefactors of the cubic and dominant quar- 
tic terms change gradually along the transition line in the 
concentration-temperature phase diagram. Either of them 
may vanish at certain points (c„, T„), where new transition 
lines originate. The new phases are classified. 
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A simple mixture of a metallic host with a magnetic 
atom, such as Aui_ c Fe c , is known to have a rather com- 
plicated phase diagram. According to Mydosh 0] the 
following phases occur at zero temperature: At very 
low c there is the Kondo-regime of independently com- 
pensated spins in a metallic host. At somewhat larger 
concentrations there is a spin glass phase of interacting 
single spins with T g oc c. For 0.5% < c < 10% the 
spin glass experiences gradual cluster formation, while 
for 10% < c < 16% one has the cluster glass phase. For 
c > 16% one enters the percolated ferromagnetic phase, 
which partly also behaves as a cluster glass. 

The Kondo regime, and the low concentration spin 
glass phase are relatively well understood. The latter is 
described by an Edwards- Anderson model with RKKY 
interactions. Its properties are obtained from a mean 
field approach |^| and from numerical analysis, see e.g. 
||. Whether or not a thermodynamic phase transition 
occurs in zero field or even in non-zero field remain topics 
of much controversy. 

Though ferromagnetism by itself is well known, clus- 
tering properties of inhomogeneous ferromagnets are also 
far from well understood. It is known that replica sym- 
metry breaking may occur before the onset of ferromag- 
netism, H possibly describing Griffiths singularities. 

The situation for the clustering spin glass (with clus- 
ters containing up to five atoms) and the cluster glass 
(where as many as 2000 atoms may build a cluster; these 
clusters order in a glassy way) is less satisfactory. Lit- 
tle is known about these phases. There seems to be no 
experimental evidence that the given names correspond 
to thermodynamic phases that are significantly different 
from the spin glass phase. Nevertheless, the existence of 
new glassy phases is the main question we wish to inves- 



tigate theoretically in this work. 

Recently one of us formulated a field theory for site- 
disordered Ising systems. With exception of the Kondo 
regime, this applies to the whole phase diagram of sys- 
tems such as mentioned above. We thus consider a sys- 
tem with translationally invariant pair couplings J(r—r') 
with a fraction c (0 < c < 1) of the lattice sites occupied 
at random. We restrict ourselves to the second order cu- 
mulant expansion. This description is Gaussian in the 
magnetization fields, and equivalent to a variational (" 
Hartree") approximation. It is quite close to the one of 
the SK-model since it involves only the order parameters 
Qa/3 (= (saSp) for small c) and their conjugates p a p. The 
replicated free energy per spin reads 

PFn = hj Wf ^ {ln(1 - c ^ q ^ w 
+ 2 H qa ? Pa p + Z 11 { l - tr * } exp x(l) ) 
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replicated spins, and tr£ Z ^ denotes the sum over ' 
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This expression is quite rich, and embodies the effect of 
clustering. Indeed, by expanding the logarithm in pow- 
ers of q a ^p one observes an effective coupling J c s{k) = 
J(k)/(1 — c(3J(k)qd), due to the presence of the diago- 
nal elements q aa = qd{c,T) < 1. If J is peaked at some 
k, Jcff(fc) will be peaked much stronger, thus exhibiting 
clustering effects. When J(k) = Jo for fcrj < k < hi, while 
vanishing elsewhere, one considers the long range, oscil- 
lating interaction J e s (r) ~ (fcrj cos k^r — k\ cos k\r) jr 2 at 
large r. In the scaling limit fej — fco ~ c ~ > 0, mean field 
becomes exact. Eq. ([!]) then has as limit the Hopfield 
model and the Sherrington-Kirkpatrick-model. |Q| 

^From eq. (Q) a Ginzburg-Landau free energy can be 
derived. Omitting the paramagnetic background, elimi- 
nating the g's and fluctuations of pd = p a a, and denoting 
Pa/3 again by q a p, we end up with 
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where now q aa = and h 2 — (3 2 H 2 /i2- The prefactor of 
the quadratic term, r = ((122 — T 2 / CJ2) /2, vanishes at the 
spin glass temperature T g (c) = V C ^2A 1 22- Furthermore, 
w = V222+T 3 J 3 /(cJ$), yi = 3^2222/2+^44/6-^422, 2/3 = 
M2222+7 l4 (>/2'/4 — 2J|)/(cJ|), and a similar expression for 
2/2- We introduced the moments of the effective coupling 
Ji = J (27r)': j [Jcs(k)) 1 and the spin-moments 
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The paramagnetic behavior is coded in the parameters 
Pd and qd, that satisfy the coupled mean field equations 
Pd = (iJ\ and qd = (J-2- All information on clustering is 
contained in r, u>, and the y's, so in the /i's and the J/. In 
the limit c — > the /i's go to unity and for T ~ the J3 
and J4 terms vanish, so that one recovers the Ginzburg- 
Landau free energy of the SK-model. The important fac- 
tors then are w = 1, 2/1 = 2/3, while the values of 2/2 
(= —2) and 2/3 (= f ) are irrelevant. When following the 
transition line T = T g {c) in the c — T phase diagram as 
function of c, it is seen that the higher /i's are rapidly os- 
cillating functions. For instance, if J 3 / J^ 2 ~Ja/ J% ~ 0, 
then yi changes sign at c = 2.7% and at c = 4.3%, while 
w becomes negative at 6.7%. 

Based on these observations we are led to assume that 
the relevant physics near the phase transition(s) is still 
contained in the GL free energy (Q). However, there is no 
reason to assume that w and y\ will always be positive. 
(A sign change of 2/1 occurs also in a Potts glass. ||) 
Given the type of the lattice and the values of the spin- 
spin couplings, the c — T phase diagram may exhibit a 
limited number of special points (c*,T*) where either w 
or y\ vanishes, and new phase boundaries originate. 

When the q a @ are expressed in the Parisi order param- 
eter function q(x), one obtains the following free energy: 
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J dx^q(x) + ^q 2 (x)-^q(x)T(x) 
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with T(x) = xq 2 (x)/2 + q(x) £ dyq(y) + dyq 2 (y)/2. 

We first investigate the region where w goes through 
zero while y\ > is fixed. In figure 1 we 

depict a fictitious phase diagram with such a situation. 
On the side where w > one has the well known spin 
glass solution of the Parisi type, as depicted in figure 2a. 
The interesting domain is w < and r ~ w 2 , since 2/2 and 
ys become relevant. In order to find an acceptable solu- 
tion we assume that j/3 < —y\ so that a = \/—y\/yz < 1- 
At h = the spin glass order parameter function 
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has plateau value q% = q(xi), determined by 

T = wqi - ~(yt + y 3 )q 2 + y (1 - x x )q\ + y/ 2 (6) 

where 1% = J* 1 dyq 2 (y). The solution is physically ac- 
ceptable as soon as y<i exceeds a certain bound and ex- 
ists for parameters such that x\ ranges from X\ = up 
to x\ = a. For (c, T) such that x± — > a the solution 
squeezes and becomes a IRSB solution with the lower 
plateau at go = (in zero field), see figure 2b. 

For w < IRSB solutions are present in a whole do- 
main. In general, a IRSB occurs in two shapes, static 
and dynamic. The static case describes physics on ex- 
ponentially large time scales where the system can over- 
come the free energy barriers between pure states. Here 
one maximizes the free energy wrt x\, which yields the 
plateau value reads 
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It sets in from x\ = 1 as a first order phase transition 
without latent heat at temperature 
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IRSB 



Tg(c) -Tg= Tg( C ) 
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Whereas the transition from paramagnet to spin glass 
has a continuous specific heat, the analogy to real glasses 
makes us expect that (also beyond mean field) the spe- 
cific heat jumps downwards at the transition PM— >1RSB. 
Both the SG and IRSB phases will exhibit a difference 
between field cooled and zero field cooled susceptibilities. 

In mean field the metastable states have infinite life- 
time. Therefore the dynamical IRSB equations lead to 
a sharp phase transition at temperature T c > T g . The 
thermodynamics of this dynamical transition is uncom- 
mon. H The entropy of the frozen state is much below 
the paramagnetic one. A crucial role is played by the 
complexity (configurational entropy), which is extensive. 
This scenario explains thermodynamically why the dy- 
namical glass transition takes place: the system just goes 
to the available state with lowest free energy. JhJ] Beyond 
mean field the dynamical aspects are reflected in the de- 
pendence on the cooling rate. 

For a dynamical IRSB-phase the qi-plateau is 
marginally stable and equal to 
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22/1 + 2/3^i(3 - xi) 
This dynamical solution sets in at a larger temperature 

^ SB -T ff (c)-r c ^r g ( C ) + ^-^ (10) 
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Both the static and dynamical solutions exist down to 
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T„ M= T,(c)-r., S T l{ c)-^-(l + ^- j ) 
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This is exactly the line where, coming from positive w, 
the SG solution gets squeezed into a 1RSB solution. The 
full phase diagram is depicted in figure 3. 

Next we consider the situation where 2/1 goes through 
zero, while w > is fixed. (In case w < the system 
will already have undergone a non-perturbative first or- 
der transition at some negative t.) One now expects 
a transition from a spin glass phase phase (y\ > 0) to 
a replica symmetric or Edwards-Anderson (EA) phase 
(yi < 0). In the EA phase there is no difference between 
field cooled and zero field cooled susceptibility. 

As it was the case for Parisi's solution of the SK- 
model, the values of 2/2 and 2/3 are now irrelevant. How- 
ever, higher order replica symmetry breaking terms will 
become relevant. All fifth order terms have been con- 
sidered for the above model. The most dangerous one 

is -(2/5/8) Y.aplaP^ 2 )^ with V5 = 6/i22222 ~ V 4 222 + 

§M442- (F° r SK: 2/5 = 8/3). We can absorb this term in 
our previous free energy using the saddle point equation 
(q 2 ) a f3 ~ —2Tq a j3/w : which amounts to replacing j/i by 
III = V\~ ^ T yblw. The most dangerous sixth order term 
is -(2/6/6)^0/3^/3, where y 6 = ^222222 - ^M42222 + 

Y§M4422 + 3^6222 ~ §M642 + 247)^66 (V6 = 16/15 for SK). 

The interesting region is where the q^p term is of same 
order of magnitude as the g^ term. This occurs when 
2/4 = jjiw 2 /2t 2 is of order unity. At fixed small positive r 
we now follow the system by changing 2/4- We thus vary 
c and T over a line at fixed distance r to the critical line. 
This is indicated by the dotted line in figure where w 
should now read 2/1 • For j/4 > 1 we will have a standard 
SG, while for 2/4 <C —1 there is the EA phase. 

When j/6 > is fixed, we find that in between the 
SG phase and the EA phase there is a SG phase with 
go > 0, although there is no external field. Coming from 
the SG-phase, go starts to become non-zero at 2/4 = 0~. 
For 2/4 — > —22/6 replica symmetry is restored since qo 
approaches q\. The y% — t phase diagram for the case 
2/6 > is shown in figure [|. As it is the case with the 
AT-line in a field, the transition EA — > SG (go 7^ 0) may 
very well be smeared beyond mean field. 

When 2/6 < we find a new, discontinuous order pa- 
rameter function, that we call SG III: q(x) — q c (x) for 
x < x\, while q{x) = q\ > q c {x\) for x > x\, see figure 
2c. As for static 1RSB solutions, the plateau has stable 
fluctuations. Coming from the EA-phase, SG III sets 
in with x\ = 0, leading to irreversibility. With respect 
to the EA-phase, the SG III phase has a smaller replica 
free energy with a discontinuous slope. There occurs a 
static first order transition without latent heat but with 
a discontinuity in the specific heat, as usual for glasses. 

At 2/4 = 1 1 2/6 1 the discontinuity of q(x) disappears and 



the standard SG solution takes over, see figure |5[ 

There are also other solutions with free energy between 
the ones of the EA and the SG III states. At 2/4 = 
|j/e| a 1RSB solution with marginal lower plateau occurs, 
as in a Potts glass. || Now the breakpoint sets in from 
x\ = 0. This 1RSB solution becomes unstable at 2/4 = 
3 1 j/6 1, where the go plateau is lifted and a foot grows 
near x = 0. We call this the SG IV solution, see figure 
2d. Like the SG III it exists up to 2/4 = 10 1 x/s 1 7 where 
the SG IV discontinuity disappears and it matches the 
standard SG solution (see Fig. 2a). In analogy with the 
marginal 1RSB solution, we anticipate that this 1RSB- 
SG TV traject is the one that occurs in dynamics. 

Also in the standard region where w and 2/1 are still 
positive some clustering effects occur. Consider the slope 
of the field-cooled susceptibility xfc = P(l — f dxq(x)). 
At T~ one has d XF c/dT = -T^ 2 + (wTg)- 1 dr / dT . In 
mean field models with 00RSB xfc is usually constant 
below T g , so these two terms cancel. There does not 
seem to be a general reason for this. Experimentally, 
the values in the SG-phase are usually lower than at T g . 
However, in the mechanically milled amorphous Co2Ge 
spin glass of Zhou and Bakker, that has about 67% of 
magnetic atoms, one expects large clustering effects. In- 
deed, xfc is monotonically decreasing with T. jll| Both 
these phenomena can be explained by our formula. 

So far our results concern mainly mean field. Whether 
or not fluctuations change them qualitatively is unknown. 

In conclusion, we have proposed a Ginzburg-Landau 
free energy for site-disordered spin glasses. It is moti- 
vated that the prefactors of the cubic and quartic terms 
can have zeroes. From these points new transition lines 
originate. We find spin glass phases of the Parisi type 
(00RSB), with 1RSB, without RSB (EA-phases), and of 
new types, the SG III and SG TV phases. For the latter 
phases the dynamics will be of new nature. 
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FIG. 1. c — T phase diagram for a fictitious system with a 
line w(c, T) = 0. PM=paramagnet; SG=spin glass. 
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FIG. 2. Shapes of the spin glass order parameter function, 
a) standard form for infinite order replica symmetry break- 
ing; b) one step replica symmetry breaking solution, c) the 
discontinuous SG ///function ; d) the SG IV function. 




FIG. 3. t — w phase diagram for a system with yi > 0, 
2/3 < — yi and j/2 sufficiently positive; with w increasing from 
right to left it may appear in Fig. 1 around the point (c, T*). 
The full (dashed) lines are static (dynamical) transition lines. 
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FIG. 4. j/i — t phase diagram for w > 0, ye > 0. The 
function q(x) in the SG phase is drawn in figure 2. a for the 
case qo = 0. In the EA-phase q(x) is constant (no RSB). 



PM 

o yi 




io|y 6 | 



FIG. 5. yi — t phase diagram for w > 0, ye < 0. In the 
SG/// phase q(x) is as in figure 2c. Dynamically this phase 
splits up in a 1RSB phase and a SG IV phase , see fig. 2b,d. 
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